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Dissipative one-dimensional maps may exhibit special points (e.g., chaos threshold) 
at which the Liapunov exponent vanishes. Consistently, the sensitivity to the initial 
conditions has a power-law time dependence, instead of the usual exponential one. The 
associated exponent can be identified with 1/(1 — q), where q characterizes the nonexten- 
sivity of a generalized entropic form currently used to extend standard, Boltzmann-Gibbs 
statistical mechanics in order to cover a variety of anomalous situations. It has been 
recently proposed [Lyra and Tsallis, Phys. Rev. Lett. 80 (1998) 53] for such maps the 
scaling law 1/(1 — q) = l/ct min — l/a max , where ct min and ct max are the extreme values 
appearing in the multifractal f(a) function. We generalize herein the usual circular map 
by considering inflexions of arbitrary power z, and verify that the scaling law holds for a 
large range of z. Since, for this family of maps, the Hausdorff dimension df equals unity 
Vz in contrast with q which does depend on z, it becomes clear that df plays no major 
role in the sensitivity to the initial conditions. 
PACS Number (s): 05.45.+b, 05.S0.-y, 05.70. Ce 
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I. INTRODUCTION 



Whenever a physical system has long-range interactions and/or long-range microscopic memory and/or 
evolves in a (mult i) fractal-like space-time, the extensive, Boltzmann-Gibbs (BG) statistics might turn out 
to be inadequate in the sense that it fails to provide finite values for relevant thermodynamical quantities 
of the system. In order to theoretically deal with nonextensive systems of this (or analogous) kind, two 
major formalisms are available: the so-called quantum groups |l[ and the generalized thermostatistics 
(GT), proposed by one of us a decade ago p|. This two formalisms present in fact deep connections 0. 
We focus here the GT. Within this framework, noncxtensivity is defined through a generalized entropic 
form, namely 

s q = k ^= lP * (q G Ti) (1) 

where k is a positive constant and {p^ is a set of probabilities associated to W microscopic configurations. 
We can immediately check that the q — > 1 limit recovers the usual, extensive, BG entropy — Yli=i Pi ^ n Pi- 
Also, if a composed system A + B has probabilities which factorize into those corresponding to the 
subsystems A and B, then S q {A + B)/k = S q (A)/k + S q {B)/k + (1 - q)S q {A)S q {B)/k 2 . This property 
exhibits the fact that q characterizes the degree of nonextensivity of the system. 

Especially during the last five years, a wealth of works have appeared within this formalism. In fact, 
it is possible to (losely) classify these works as follows: (i) Some of them j| address the generaliza- 
tion of relevant concepts and properties of standard thermostatistics, such as Boltzmann's H-theorem, 
fluctuation-dissipation theorem, Onsager reciprocity theorem, among others; (ii) Other GT works || 
focuse applications to some physical systems where BG statistics is known to fail (stellar polytropes, 
turbulence in electron-plasma, solar neutrino problem, peculiar velocities of spiral galaxies, Levy anoma- 
lous diffusion, among others), and yields satisfactory results; (iii) Finally, an area of interest which is 
progressing rapidly, addresses the long standing puzzle of better understanding the physical meaning of 
the entropic index q. This line concerns the study of nonlinear dynamical systems (both low ^-^] and 
high ]lO| , pT| dimensional dissipative ones, as well as Hamiltonian systems O]) in order to clarify the 
connection between q, the sensitivity to the initial conditions and a possible (multi)fractality hidden in 
the dynamics of the system. This paper belongs to the last class of efforts and is organized as follows. 
In Section 2 we briefly summarize recent related results. In Section 3 we introduce a new map which 
generalizes the circular one, and study its main properties. Finally, we conclude in Section 4. 



II. POWER-LAW SENSITIVITY TO INITIAL CONDITIONS 



The most important dynamical quantities that are used to characterize the chaotic systems are the 
Liapunov exponent Ai and the Kolmogorov-Sinai entropy K\ (the meaning of the subindex 1 will soon 
become transparent). Let us define, for a one-dimensional map of the real variable x, the quantity 
£(t) = limAx(o)^o Aa(o) ' wnere Ax(0) and Ax(t) are discrepancies of the initial conditions at times 
and t respectively. It can be shown that, under quite generic conditions, £ satisfies the differential 
equation d^/dt = Ai £, hence = exp(Ait). Consequently, if Ai < (Ai > 0) the system is said 
to be strongly insensitive (sensitive) to the initial conditions. Similarly, for a dynamical system under 
certain conditions, we can define K\ as essentially the increase, per unit time, of S\ = — Y] : _ 1 Pi In p^. 
Furthermore, it can be shown that, with some restrictions, K\ = Ai if Ai > and K\ — otherwise. 
This is frequently referred to as the Pesin equality |l3| . 

The case we are focusing in the present work is the so called marginal case, corresponding to Ai = 0. 
It has been argued that, in this case, the differential equation satisfied by £ is d^/dt = X q £ 9 , hence 

£(i) = [l + (l-9)V] V(1 " 9) . (2) 

where X q is the generalized Liapunov exponent. One can verify that q = 1 recovers the standard, 
exponential case whereas q ^ 1 yields a power-law behavior. If q > I (q < 1) the system is said to 
be weakly insensitive (sensitive) to the initial conditions. Furthermore, for this marginal case we can 
define the generalized Kolmogorov-Sinai entropy K q as the increase, per unit time, of S q . Finally, for 
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this anomalous case, it has been argued M that the Pesin equality itself can be generalized as follows: 
K q = \ q if X q > and K q = otherwise. 

Recently, these ideas have been applied to some dissipative one-dimensional maps (a logistic-like and 
a periodic-like map, sharing the same universality class, as well as the standard circular map, which 
belongs to a different universality class), and the numerical results suggested a close relationship between 
the nonextensivity parameter q and the fractal (Hausdorff) dimension df associated with the dynamical 
attractor 0-^]. Very specifically, in those examples, when df approaches unity (which is the Euclidean 
dimension of the system) from below then q also approaches unity from below, and does that in a 
monotonic manner. Naturally, this fact strongly suggests that the validity of the statistical q — 1 (BG) 
picture is intimately related to the full occupancy of the phase space. However, an important question 
which remains open is whether the full occupancy is sufficient for having a BG scenario. We will show 
here that it is not! We shall introduce in this paper a generalized circular map (characterized by an 
inflexion power z), and show that, at the critical point, df(z) = 1 (Vz), and nevertheless q < 1 ! This 
fact consistently complements the new scaling relation proposed in M , namely 



t q (Xmin &max 

where a m j„ and a max are the extremes of the multifractal singularity spectrum /(a) of the attractor (for 
details see |ljl). This relation clearly indicates that once the scaling properties of the dynamical attractor 
are known, one can precisely infer the proper entropic index q that must be used for other purposes. 



III. A FAMILY OF CIRCULAR-LIKE MAPS 



The circle map is an iterative mapping of one point on a circle to another of the same circle. This map 
describes dynamical systems possessing a natural frequency u)\ which are driven by an external force of 
frequency w 2 ; ^ = uji/u>2 is known as the "bare" winding number. These systems tend to mode-lock at 
a frequency uj* and u) = uj\/u)2 is known as the "dressed" winding number. The standard circle map, for 
one-dimension, is given by 

0*+i = fi + t - — sin(27T0 t ) mod(l), (4) 

with 0<f2<l;0<-K"<oo. For K < 1 the circle map is linear at the vicinity of its extremal point and 
exhibits only periodic motion. From now on we take K = 1, the onset value above which chaotic orbits 
exist. For this map, once mode-locked, u) — ]irat— >00 (dt+i ~ &t) remains constant and rational for a small 
range of the parameter £1 with the "dressed" versus "bare" winding number curve exhibiting a "devil 
staircase" aspect (if t+1 < 6 t then one shall use u = 1 + 6 t+ i — 9 t in order to leave it mod(l)). At 
the onset to chaos, a set of zero measure and universal scaling dynamics is produced at special irrational 
dressed winding numbers which have the form of an infinite continued-fraction expansion 

u> = — ■ — — , (5) 



with n,m,p, .. integers. The best studied one is when u> equals the golden mean, i.e., ujqm — (v5 — l)/2 
(n = m = p = ... = 1) [ |il^ , |l7[ . It is worth mentioning that the golden mean is the asymptotic ratio 
between consecutive numbers of the Fibonacci series (limn^oo F n /F n+ i = (y5 — l)/2, where Fq = 0, 
Fi = 1 and F n — F n -\ + F n _2). In order to determine the bare winding number at this critical 
point, we iterate the map for a large number of time steps (10 6 steps) starting with 8q = and use a 
linear regression to numericaly compute to. The bare winding number 11 is then adjusted to have the 
renormalized winding number u> equal to the golden mean, which results in tt c = 0.606661... With these 
parameters, the standard circle map has a cubic inflexion {z = 3) near its extremal point 8 = 0. 
A generalized version of the circle map can be defined as 
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where z > (z — 3 reproduces the standard case). For every value of z, the golden mean of to corresponds 
to different "bare" winding numbers, which we call as f2 c . In order to determine these critical values of 
D,, one searches, within a given precision (12 digits in our calculations), the value of fl corresponding to 
logm with the same precision. The calculated values of Q c , as a function of z, are shown in Fig. 1. The 
numerical values are indicated in the Table. We remark that, in the limit z — > [z — > oo) we verify that 
n c oc z 1 / 2 (1 - C oc 1/z 13 with (3 ~ 0.41). 

For our present purpose, a very important feature of this map is that, for every value of z, the critical 
attractor visits the entire circle (0 < 9 t (mod 1) < 1) and therefore has a support fractal dimension 
df = 1. The calculated values of df, using a box counting algorithm, are indicated in Fig. 2. In what 
concerns the sensitivity to the initial conditions, the function £(t) is given by 



lnf(t) =ln 




— sin(27T0 t ) 

Z7T 
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t - -!- sin(27T0 t ) 



[1 - COs(27T0 t )] } (7) 



and displays, for f2 = f2 c , a power-law divergence, £ oc ^/t 1 -?) , from where the value of q can be calculated 
by measuring, on a log-log plot, the upper bound slope, 1/(1 — q). In Fig. 3, the z — 4.5 and z = 6 cases 
have been illustrated; see also the Table. 

In order to check the accuracy of the scaling relation (3), one needs to determine the a m in and a max 
values of the f(a) curve. Therefore, one has to study the structure of the trajectory 6\, 02, and 
to estimate the singularity spectrum (strength of singularities a and their fractal dimensions /) of this 
Cantor-like set. To perform the numerical calculation we truncate the series 9i at a chosen Fibonacci 
number F n (we recall that F n /F n+ i gives the golden mean and therefore defines the proper scaling 
factor). The distances U between consecutive points of the set define the natural scales for the partition 
with measures pi — l/F n attributed to each segment. After that, the singularity spectrum can be 
directly obtained following a standard prescription jL4|. In general, sup Q f(a) = df and, in the present 
case, f{a min ) = f(a max ) = 0. We define a top through f(a top ) = df. 

However, the situation for the generalized circle map is somewhat different than that of the standard 
circle map in the sense that the standard one has a fast convergence of the f(a) curve when larger 
number of iterations are considered, whereas the generalized map presents only a slow and oscillatory 
convergence. An example of a sequence of f(a) curves obtained from increasing number of iterations is 
shown in Figure 4. Notice its non-monotonic behavior, specially near the upper edge. In Figure 5, we plot 
the numerically obtained values of a m i n , atop and Oi max as a function of 1/ IniV, where N is the number 
of iterations. From these data, we are not able to accuratelly estimate their asymptotic values for large 
map inflexion z. In Figure 6 we plot the extrapolated f(a) curves for typical values of z. Although these 
show the main expected trends of the singularity spectra, namely z-dependent shape but df = 1 for all z, 
their extremal points may need further corrections. We shall point out that this feature is inherent to the 
numerical method used to estimate the /(a) curve representing the singularity strengths of a multifractal 
measure. Its extremal points are governed by the scaling behavior of the most concentrated (a m i„) and 
most rarefied (a max ) sets in the measure, the latter being usually poorly sampled. 

An alternative method for computing the extremal values of the singularity strengths a can be obtained 
by studying how the distances around 9 = scale down as the trajectory 9i is truncated at two consecutive 
Fibonacci numbers, F ni F n+ \. Shenker has found that this distance shall scale by a universal factor cxf(z) 
|l6f (F stands for Feigenbaum). This region around 9 = corresponds to the most rarefied one so that 
Z_oo ~ [af(z)] _ ™. The corresponding measure scales as p~oo = Pi = 1/F n ~ (ugm)", which leads to 

§00 

lnp^op _ In lo gm fR , 

ln/_oo ln[OF(z)J 

Following along the same lines, the most concentrated region on the set shall scale down as Z+oo ~ 
ap(z)~ zn while p +OQ ~ [wgm]", so that 

__ lnp +co In logm /„. 

Oirain — . . — , , , • \s)) 

m/ +00 Ya.\a F \z)\ 2 

Eqs. (8) and (9) imply 
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(10) 



In Figure 7, we show the critical sequence of the distances around 6 — for a large value of the map 
inflexion (z = 8). The data provide an accurate estimation of the universal factor oip (z — 8) = 1.1568 
from which precise values of a rn i n and a max can be inferred. In the Table we list the results for 3 < z < 8, 
together with the values of a>t op from the extrapolated f(a) curves and q obtained from both the scaling 
relation (3) and from the sensitivity function (7). Notice that the scaling relation (3) is satisfied for all 
z (see also Figure 8), just like the case of the logistic-like maps ||. Furthermore, these results indicate 
two other important points: (i) Even though the fractal dimension of the support df is Euclidean (i.e., 
dt = 1) for all z, the system sensitivity to initial conditions is still z-dependent. (ii) What matters is 
a m i n and a max , and not df, in other words, what precisely controls the entropic index q is not df but 
the sensitivity to initial conditions, reflected by the multifractal nature of the attractor. 



IV. CONCLUSIONS 



In this paper, we contribute to the field of low-dimensional dissipative systems by introducing a conve- 
nient generalization (with inflexion power z) of the standard circular map and then studying its critical 
point (analogous to the chaos threshold of logistic- like maps). More precisely, we have numerically stud- 
ied its sensitivity to the initial conditions and have shown that it is given by a power-law (£ oc t 1 '^ 9 ') 
instead of the usual exponential behavior. Moreover, we have shown that this example, as the logistic-like 
maps, satisfies the scaling law given in Eq. (3). Although q, a m i n and a max depend on z, df does 
not. This is a quite important result because it illustrates that, for having a Boltzmann-Gibbs scenario 
(q = 1), it is not enough to fully occupy the phase space during the dynamical evolution of the system. 
What is essentially necessary is to have a quick, exponential-like occupation of the phase space, so that 
ergodicity and mixing are naturally attained. 

In addition to this, it is worth mentioning that there are also other efforts along this line which address 
high-dimensional dissipative systems, namely those exhibiting self-organized criticality |19|. Amongst 
them, the study of the Bak-Sneppen model for biological evolution Jh]] and the Suzuki-Kancko model for 
the battle of birds defending their territories |ll[] can be enumerated. In both cases it is shown that, at 
the self-organized critical state, a power-law sensitivity to initial conditions emerges, like in the present 
case. 

As a final remark, it should be emphasized that these ideas seem to be valid and applicable not only 
to low- and high-dimensional dissipative systems but also to conservative (Hamiltonian) systems with 
long-range interactions . This fact has been illustrated very recently jl2| on the long-range classical XY 
ferromagnetic model, whose entire Liapunov spectrum collapses (for an infinitely wide energy interval) to 
zero at the thermodynamic limit if (and only if) the range of the interactions is sufficiently long. Further 
efforts focusing, along these lines, both dissipative and conservative systems are welcome. 
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Table and Figure Captions 



Tabic : Our best numerical values for il c , ap, a m i n , oitop, oi max and q from both the scaling relation 
(3) and from the sensitivity function (7) for various z. It is worthy to mention that we numerically verify 
that, for z > 3, a top > \a m in flmai] 1 ' 2 (the equality appears to hold for z = 3). 



Figure 1 : The values of the "bare" winding number U c as a function of z for the generalized circle 
maps. 



Figure 2 : Box counting graph for determining the fractal dimension dj of the critical attractor for 
typical values of z. 



Figure 3 : The plot of ln£(7V) versus In TV. (a) for z = 4.5 and (b) for z = 6. 



Figure 4 : Approximate multifractal singularity spectra of the critical attractor of the z = 8 generalized 
circular map for increasingly large number of iterations, going from N = 233 to 987. 



Figure 5 : a min (bottom three lines), a top (middle three lines) and a max (top three lines) from the 
singularity spectra obtained from distinct number of iteractions and typical values of z. Notice the slow 
and oscillatory convergence. 



Figure 6 : The extrapolated f(a) curves for typical values of the map inflexion z. Notice that, 
although the shape is z-dependent, they present df = 1 for all z. The solid lines correspond to f(a) = 1 
and f(a) = a. The dotted lines are guides to the eye. 



Figure 7 : The sequence Ot as a function of t for z = 8. The minimal distance to 8 = scales down as 
£-0.305 Tj s i n g that F n <~ for large n one obtains that Z_oo <~ ap n , with ap(z = 8) = 1.158. 



Figure 8 : l/a min — l/a max versus 1/(1 — q) values from the sensitivity function. The straight line 
represents the scaling prediction (Eq. (3)). 
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Table 



z 


n c 






Ottop 




q 

(Eq.(3)) 


q 

(Eq.(7)) 


3.0 


0.606661063469... 


1.289 


0.632 


1.096 


1.895 


0.05 ±0.01 


0.05 ±0.01 


3.5 


0.629593799039... 


1.258 


0.599 


1.124 


2.097 


0.16 ±0.01 


0.15 ±0.01 


4.0 


0.648669091983... 


1.234 


0.572 


1.167 


2.289 


0.24 ±0.01 


0.24 ±0.01 


4.5 


0.664861001064... 


1.218 


0.542 


1.213 


2.440 


0.30 ±0.01 


0.30 ±0.01 


5.0 


0.678831756505... 


1.205 


0.516 


1.266 


2.581 


0.36 ±0.01 


0.36 ±0.01 


5.5 


0.691048981515... 


1.195 


0.491 


1.314 


2.701 


0.40 ±0.01 


0.40 ±0.01 


6.0 


0.701853340894... 


1.185 


0.473 


1.351 


2.838 


0.43 ±0.01 


0.44 ±0.01 


7.0 


0.720182442561... 


1.170 


0.438 


1.451 


3.065 


0.49 ±0.01 


0.50 ±0.01 


8.0 


0.735233625356... 


1.158 


0.410 


1.518 


3.280 


0.53 ±0.01 


0.53 ±0.01 
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